Universality in the off-equilibrium critical dynamics of the 3c? diluted Ising model 
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We study the off-equilibrium critical dynamics of the three dimensional diluted Ising model. We 
compute the dynamical critical exponent z and we show that it is independent of the dilution only 
when we take into account the scaling-corrections to the dynamics. Finally we will compare our 
results with the experimental data. 

PACS numbers: 05.50.+q, 75.10.Nr, 75.40.Mg 

The issue of Universality in disordered systems is a 
controversial and interesting subject. 

Very often in the past it has been argued that critical 
exponents change with the strength of the disorder [Q . 
While, on a deeper analysis, it has turned out that those 
exponents were "effective" ones, i.e. they are affected by 
strong scaling corrections. So, when one studies the crit- 
ical behavior of a disordered system it is mandatory to 
control the leading correction-to-scaling in order to avoid 
these effects that could modify the dilution-independent 
values of the critical exponents. For instance, in Rcf. ||] 
the equilibrium critical behavior of the three dimensional 
diluted Ising model was studied. The authors showed 
that taking into account the corrections-to-scaling it was 
possible to show that the static critical exponents (e.g. v 
and rf) and cumulants were dilution-independent. These 
numerical facts supports the (static) perturbative renor- 
malization group picture: all the points of the critical line 
(with p < 1) belong to the same Universality class (with 
critical exponents given by the random fixed point) ||. 
Their final values of the exponents Q were in very good 
agreement with the experimental figures (see below). 

We will show that an analogous effect also happens 
in the off-equilibrium dynamics of the diluted ferromag- 
netic model and we will take it into account in our data 
analysis, in order to get the best estimate of the critical 
dynamical exponent. 

The critical dynamics of the diluted Ising model 
has been studied experimentally in Ref. H using 
neutron spin-echo inelastic scattering on samples of 
Feo.46Zno.54F2 (antiferromagnetic diluted model) and has 
been compared with the results obtained in pure sam- 
ples (FeF2) j|. For the pure model a dynamical criti- 
cal exponent z = 2.1(1) was found (in good agreement 
with the theoretical predictions based on one-loop per- 
turbative renormalization group (PRG) §) whereas m 
the diluted case the exponent z = 1.7(2) was computed 



(three standard deviations away of the analytical predic- 
tion based on (one-loop) PRG that provides z ~ 2.34 ||). 
Furthermore, the dynamical exponent was computed in 
the framework of the PRG up to two loops and it was 
obtained z = 2.237 @ and z = 2.180 f§ (the experimen- 
tal value is at 2.5 standard deviation of the two loops 
analytical result). 

In the experiment Q were measured critical ampli- 
tudes 100 times smaller than those computed in the pure 
case. It is clear that a more precise experiment on this 
issue will be welcome. We should point out that the crit- 
ical dynamics of a diluted antiferromagnet is the same as 
of a diluted ferromagnet. 

A numerical study of the on-equilibrium dynamics in 
diluted systems was performed in 1993 by Heuer |J. He 
measured the equilibrium autocorrelation functions for 
different concentrations and lattice sizes. The autocorre- 
lation time (t) depends on the lattice size (L) via the for- 
mula t oc L z (neglecting scaling corrections). He found 
that all the data, for concentrations not too close to 1, 
were compatible, for large L, with the assumption of a 
single dynamical exponent, different from the one of the 
pure fixed point and similar to the analytical estimate of 
Ref. J(| (z ~ 2.3). The final value reported by Heuer was 
z = 2.4(1). 

The main goal of this work is to check Universality 
in the critical dynamics of diluted models (i.e. whether 
the dynamical critical exponent is dilution independent) 
in the off-equilibrium regime fiofl . To do this we mon- 
itor scaling corrections in the same way it was done in 
the static simulations Therefore we will also obtain 
the value of the corrections-to-scaling exponent for the 
dynamics. Our motivation to study the off-equilibrium 
dynamics instead of the equilibrium one is based on two 
reasons. The more important reason is that the exper- 
imental data was obtained in the off-equilibrium regime 
and the second one is that (in general) it is easier to sim- 
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ulate systems in the off-equilibrium regime. Moreover it 
will be possible to confront our z computed in the off- 
equilibrium regime with that obtained at equilibrium || . 

The relevance of the corrections-to-scaling is twofold. 
The first one is that the scaling-corrections are very im- 
portant in the right determination of the static (equilib- 
rium simulation) critical exponents [Q. In some models 
the corrections-to-scaling change the anomalous dimen- 
sion of the order of 10 % (see for example Ref. The 
second one is that the correction-to-scaling exponent can 
be (and it has been) computed in a real experiment [ jl2| . 

We have studied the three dimensional diluted Ising 
model defined on a cubic lattice of size L and with Hamil- 
tonian 



H = - ^ e i e 3 S i S 3 
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where Si are Ising spin variables, < ij > denotes sum 
over all the nearest-neighbor pairs and ti are uncorre- 
cted quenched variables which are 1 with probability p 
and zero otherwise. 

We have measured, at the infinite volume critical point 
and for several concentrations p, the non-connected sus- 
ceptibility, defined by 



— ^2 ( S i S j 
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where the brackets stand for the average over different 
thermal histories or initial configurations and the hori- 
zontal bar for an average over the disorder realizations. 
The indices i and j run over all the points of the cubic 
lattice. In practice we use a large number of disorder re- 
alizations (Ns — 512) each with a single thermal history, 
what amounts to neglect the angular brackets in Eq.(||). 
This procedure is safe and does not introduce any bias. 

With the notation of the book of Ma || we can write, 
for instance, the following equation for the response func- 
tion, under a transformation of the dynamical renormal- 
ization group (RG) with step s, 
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where uj is the frequency, k is the wavelength vector, z is 
the dynamical critical exponent, by fi we denote all the 
parameters of the Hamiltonian, fj,* is the fixed point of 
the renormalization group transformation, £ is the static 
correlation length and finally y\ is the relevant eigenvalue 
(equals to \jv: y\ is the scaling exponent associated with 
the reduced temperature), e\ is its associate eigenvector 
and y 2 is the greatest irrelevant eigenvalue (y 2 < 0) of the 
renormalization group transformation (we have assumed 
that the system posses only one relevant operator). 

Using Eq.(|3|) and considering the leading scaling- 
corrections for a very large system at the critical 



temperature, we can write the dependence of the suscep- 
tibility on the Monte Carlo time as 



(4) 



where t is the Monte Carlo time, T c (p) is the critical tem- 
perature, A(p) and B{p) are functions that depend only 
on the spin concentration, 7 is the exponent of the static 
susceptibility, v is the exponent of the static correlation 
length, z is the dynamical critical exponent and finally 
w = — 2/2 is the correction-to-scaling exponent. Hereafter 
we denote W4 = w/z. We recall that w corresponds with 
the biggest irrelevant eigenvalue of the RG in the dy- 
namics, in principle w will be different from the leading 
correction in the static (that we will denote by w s ) 0. In 
addition, an analytical correction-to-scaling comes from 
the non singular part of the free energy and gives us a 
background to add to Eq.(|4|). In our numerical simula- 
tions we can neglect this background term (i.e. we will 
show that ^/(vz) — uj/ z ~ 0.5 » 0). 

Moreover, Eq.(^) is valid for times larger than a given 
"microscopic" time and for times (in a finite lattice) less 
than the equilibration time (that is finite in a finite lat- 
tice). 

To study numerically the present issue we have simu- 
lated L = 100 systems for different spin concentrations 
p = 1.0,0.9,0.8,0.65,0.6,0.5 and 0.4 at the critical tem- 
peratures reported in Ref. [2|. The Metropolis dynam- 
ics provides our local dynamics. We have checked 
in all the simulations that we were in an off-equilibrium 
situation: for the volumes and times we have used, the 
non-connected susceptibility is far from reaching its equi- 
librium plateau (in a finite system). For completeness we 
also report the numerical estimate of the critical expo- 
nents for the random fixed point, where all the systems 
with p < 1 should converge for large length scales g] : 7 = 
1.34(iyj/ = 0.6837(53), 77 = 0.0374(45) and w s = 0.37(6) 
(Ref. flqj , using PRG, provides to = 0.372(5) in the mas- 
sive scheme and uj — 0.39(4) in the minimal subtraction 
one). It is worth noting that experimentally the best es- 
timate of the susceptibility exponent is 7 = 1.33(2) |]l6| . 

At this point we can recall the one-loop prediction of 

the PRG for the v and r\ exponents: v = -| + ■j-^/ff 
and 77 = —e/106 PJI^|, where e = 4 — d, d being the 
dimensionality of the space. If we substitute e = 1 we 
obtain the following ("naive") estimates: 77 = —0.0094 
and v = 0.5841. Obviously the previous naive estimates 
are far from the numerical and experimental values of the 
critical exponents. This would also imply that even the 
one-loop PRG estimate of the dynamical critical expo- 
nents will stay far from the true value. 

Notice also that the anomalous dimension exponent 
(77) takes near the same value either at the pure and 
at the random fixed point. One can argue that this 
holds using the arguments provide in Ref. p8| using a 
e'-expansion (where d = 2 + e') |19|. This fact, assuming 
the naive dynamical theory (Van Hove theory or con- 
ventional theory) Q , implies that the dynamical critical 
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exponent z — 7/1/ = 2 — 7/ is the same for both, diluted 
and pure system, to first order in e'. We will show that 
this is not the case for our diluted model. The Van Hove 
theory was used in Q to interpret the experimental data. 

An analytical estimate of the value of the dynamical 
critical exponent has been taken from Ref. || where a 
dynamical -^e-expansion (e = 4 — d) was done: z = 
2 + y / 6e/53 + 0(e), that in three dimension becomes 
z ~ 2.34, where we have neglected the terms 0(e). We 
can recall the two loops computation: z — 2.237 M and 
z = 2.180 (|). One of the results of this work should be 
about the reliability of the previous estimates of z (the 
first and second term of an -^/e-expansion) . 

With all these ingredients we can analyze our numer- 
ical data for the dynamical non-connected susceptibil- 
ity and check whether or not the Universality, based on 
renormalization group arguments, holds. 
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FIG. 1. The growth of the out-of-equilibrium susceptibility 
with the Monte Carlo time, at the critical temperature. The 
lattice volume is always 100 3 and the spin concentrations are 
reported in the plot. The errors are smaller than the symbols. 



In the first plot (Fig. [j]) we show the numerical data 
in a double logarithm scale. The slope gives, neglect- 
ing the corrections-to-scaling, the ratio ^ /(vz). It seems 
that all the lines behave in a power law but with different 
slopes [i.e. different exponents j/(uz)]. This fact could 
call for non universality in this model (i.e. critical expo- 
nents vary along the critical line). In addition, if we take 
into account the main result from the static [|| which 
states that the static critical exponents (e.g. v and r\) do 
not depend on the dilution degree, we obtain a dynami- 
cal exponent that depends on the dilution, violating the 
prediction of the dynamical perturbative renormalization 
group ||. In fact, following the RG flow (for p < 1) we 
should always end at the same random fixed point and 
so, for large scales (in time and space) z is not expected 
to depend on the dilution degree. 

In the previous analysis we have not taken into account 
the scaling corrections. However, we are able to monitor 
the leading scaling-corrections given by the exponent w. 
We succeeded in fitting (using the MINUIT routine |oJ) 



all our numerical data to Eq. (|4j) for 0.5 < p < 0.8. We 
have 10 parameters to fit : A(p) and B(p) for four dilu- 
tions (p = 0.8, 0.65, 0.6, 0.5), ^j(vz) and ui/z, these last 
exponents assumed dilution independent. 

In this way we have computed the functions A(p) and 
B(p) in Eq.(Q) and ^/(vz) and uj/z. By fitting the data 
using t > 4 we have obtained a very good fit (with 
X 2 /d.o.f = 33.8/34, where d.o.f stands for degrees of free- 
dom) and the following values for the dynamical critical 
exponent and the leading dynamical scaling-corrections 



z = 2.62(7) , u = 0.50(13), 



(5) 



where we have used the value of the static critical expo- 
nents 7 = 1.34(1) and v = 0.6837(53) @. 

In order to check the stability of the previous fit we 
have tried a new fit using only times t > 8. The fit again 
is very good (with x 2 /d.o.f = 29.7/30) and 



z = 2.58(7) , w = 0.72(16), 



(6) 



Clearly the fit is very stable since both exponents are 
compatible inside the error bars (one half standard devi- 
ation in z and one standard deviation in u>). Therefore we 
take, as our final values, z = 2.62(7) and to = 0.50(13). 

In Fig. H we show our results for the amplitudes A(p) 
and B(p) (using the results of the fit with t > 4; t = 4 
plays the role of the microscopic time for this model and 
algorithm, see the previous discussion). The main re- 
sult of these fits is that the numerical data can be well 
described using dilution independent exponent (both dy- 
namical and static), while the value of the dilution only 
enters in the non-universal amplitudes, A(p) and B(p). 
This fact clearly supports Universality in this model. 
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FIG. 2. The amplitudes defined by Eq.(^) are smooth func- 
tions of the spin concentration. Where the B(p) crosses the 
axis a "perfect Hamiltonian" can be defined (see text). 

From Fig. || we can compute the value of the dilution 
in which there is not (leading) scaling-corrections (one 
kind of "perfect Hamiltonian" for this dynamical prob- 
lem). For p ~ 0.63 we obtain B(p) ~ and so with this 
dilution it is possible to measure dynamical critical expo- 
nents [e.g. ~f/(zi>) from the growth of the susceptibility, 
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(d — \jv)j{zv) from the relaxation of the energy, etc.] 
neglecting the underlying (leading) scaling-corrections. 
This dilution could be a good starting point in order to 
monitor the sub-leading scaling-corrections. 

Systems with spin concentrations p = 0.9 have also 
been simulated, but the data from these runs are not 
been included in the previous analysis, because they can 
not be well fitted with the formula of Eq.([|). We can ex- 
plain this fact assuming that for this dilution the system 
is in the cross-over region, for the lattice and times we 
used. Also in the static studies a similar effect was found 
and only for p < 0.8 was possible to obtain final values 
(for exponents and cumulants) dilution independent 0. 

In order to convince the reader of the goodness of 
our fits we plot in Fig. || the non-connected suscep- 
tibility divided by just the correction-to-scaling factor 
[A(p) + B(p) t~ Wd ]. If universality holds (i.e. all the crit- 
ical exponents, dynamical and static, are dilution inde- 
pendent) all the data points (corresponding to four di- 
lution degrees) should collapse on a straight line in a 
double logarithm scale. It is clear from this figure that it 
is what happens. The equation of the curve is t~ with 



-2- = 0.748. 
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FIG. 3. The universal part of the susceptibility growth. 
The collapse of the data for different concentrations is the 
confirmation that Universality holds. 

We have shown that it is possible to describe the off- 
equilibrium numerical data assuming critical exponents 
(dynamical as well as static) independent on the dilution 
for a wide range of dilutions. This supports the predic- 
tions of the (perturbative) renormalization group for the 
statics as well as for the dynamics. So, the (perturba- 
tive) RG scenario that predicts that all the points on the 
critical line (for p < 1) belong to the same Universality 
class is very well supported by numerical simulations. 

We have found that our estimate of the dynamical crit- 
ical exponent z = 2.62(7) is incompatible with the exper- 
imental value z = 1-7(2). Further numerical and exper- 
imental studies should be done in order to clarify this 
discrepancy. 

We can compare the value of the dynamical critical ex- 



ponent computed off- and on-equilibrium. The Hcucr's 
estimate was z — 2.4(1) and the difference with our es- 



timate z — 2.62(7) is z ff-c 



0.22(12), i.e. 1.8 



standard deviations. The conclusion is that both estima- 
tions are compatibles in the error bars. In any case, it 
will be interesting to compute z on-cquilibrium by con- 
trolling the scaling corrections. 

Moreover, our estimate is not compatible with that of 
PRG to order y/e in the -^/e-expansion (z = 2.34). The 
comparison with the two loops estimates of z jj],H is still 
worse. One possible explanation for this disagreement 
could be the lack of Borel summability the diluted model 
shows 1 21 . We remark again that the one loop PRG es- 



timates of the static critical exponents was very bad (see 
below) . 

Another interesting issue is to compare the dynam- 
ical scaling-corrections and the static ones. Unfortu- 
nately our statistical precision is unable to solve this is- 
sue. For instance, taking the values of the t > 4 we obtain 
lu — u>d — 0.13(6) that is compatible with zero assuming 
two standard deviations. If we take the values of the t > 8 
fit we obtain uj — Ud = 0.35(17). We will devote further 
work (analytical and numerical) in order to discern if the 
leading dynamical scaling-correction corresponds to the 
leading static scaling-correction. 

We wish to thank H. G. Ballesteros, D. Belanger, L. 
A. Fernandez, Yu. Holovatch, V. Martin Mayor and A. 
Munoz Sudupe for interesting discussions. 
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